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Abstract 

We present preliminary results of a search for the decays B + — > p + 7, B — > p°7, and B° — > CJ7. The 
analysis is based on data containing 347 million BB events recorded with the BABAR detector at 
the PEP-II B factory. We measure branching fractions of B{B+ -> p+7) = (1.06lo;|i±0.09) x 10~ 6 
and B(B° -> p°j) = (0.77lg;?g ± 0.07) x 10" 6 , where the first errors are statistical and the second 
systematic, and set a 90% C.L. upper limit of B(B° — > CJ7) < 0.84 x 10~ 6 . Assuming isospin 
relations between the three branching fractions, these results are used to determine the CKM 
matrix element ratio \V td /V ts \ = 0.17lig;^f (exp.^oi^theor.). 
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1 INTRODUCTION 



In the Standard Model (SM) description of the decays 1 B + — > /? + 7, B° — > p°7, and B° — > ury, the 
dominant contributions arise from b ^ dj penguin diagrams of the type shown in Figure ^ 




u, c, t 



Figure 1: Feynman diagram for b — > dq. 



Relating the three individual decay rates by isospin symmetry and using the measured ratio between 
the charged and neutral B-meson lifetimes t b +/t b q = 1.071 ±0.009 [15], one can define a combined 
branching fraction 

B[B -► (j>/w)i\ = B{B+ -► p+7) = 2^B(B° A) = 2^B(B° -► w 7 ) (1) 

The results of recent calculations of S[B — > (p/cj) 7 ] are in the range of (0.9-1.8) x 10 -6 [1-3]; 
however, these could be modified by processes beyond the SM [4]. Within the SM, the isospin 
violation in these decays is expected to be small; a recent estimate [2] is (1.1 ± 3.9)%. 

While the exclusive rates have a large uncertainty due to non-perturbative long-distance QCD 
effects, much of this uncertainty cancels in the ratio of B — > (p/co)^ and B — * K*^y rates. Since 
the dominant diagram in Figure ^ involves a virtual top quark, this ratio is related to the ratio of 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements |Vtd/Vt s | [1,2] via 



B(B -+ (p/q,) 7 ) 
B(B -» K*j) 



V td 



V, 



f.s 



• 1 - m 2 p / Ml 
1-m^/Ml 



C 2 [1 + AJ2]. (2) 



Here, the form factor ratio £ describes the flavor-SU(3) symmetry breaking between p/co and K*, 
and AR accounts for annihilation diagrams. Physics beyond the Standard Model could affect these 
decays and create inconsistencies between the measurements of |Vtd/Vt s | obtained from this analysis 
and those obtained in studies of B® and B^ mixing. 

Previous searches by BABAR [5] and CLEO [6] have found no evidence for B — > (p/uS) r ) decays. 
An observation of the decay B — > p 0/ y was recently reported by the Belle collaboration [7j. 

In this paper we report a search for the decays B + — ► /9 + 7, B° — > p 0r y, and B° — > wy. The results 
presented in this paper use a BABAR data sample containing 347 million BB events, corresponding 
to an integrated luminosity of 316 fb -1 , and supersede those in Ref. [5]. 

1 Charge conjugate modes are implicitly included throughout. 
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2 THE BABAR DETECTOR 



The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e + e~ storage ring. Charged particle trajectories are measured by a combination of a five- 
layer silicon vertex tracker and a 40-layer drift chamber in a 1.5-T magnetic field. Photons and 
electrons are detected in a CsI(Tl) crystal electromagnetic calorimeter (EMC) with photon energy 
resolution a E /E = 0.023( J E/GeV)~ 1/4 © 0.019. A ring-imaging Cherenkov detector (DIRC) is 
used for charged-particle identification. In order to identify muons, the magnetic flux return is 
instrumented with resistive plate chambers and limited streamer tubes. A detailed description of 
the detector can be found elsewhere [8]. 



3 EVENT RECONSTRUCTION AND SELECTION 

The decays B + -► p+7, B° -► p 7 , and B° — > uj are reconstructed by combining a high-energy 

photon with a vector meson exclusively reconstructed in the decays p° — > 7r + 7r~ (B « 100%), 
p+ ^+^0 ( B w 100% ) 5 and w ^ ^+^-^0 ( B = [ g9 j ± 7 ] % [ 15 ])_ 

The primary source of background is due to continuum events (e + e~ — ► gg, with q = u,d, s,c) 
that contain a high-energy photon from 7r° or 77 decays or from initial-state radiation (ISR). Decays 
of -B — > K*7, K* — ► K7r can enter the signal selection, e.g., when a is misidentified as a tt ± . 
B -> (/9/w)tt and B -> (p/a;)r? processes are also found to be relevant when a high-energy photon is 
produced in the 7r° or r\ decay. In addition, there is combinatorial background from high-multiplicity 
b — ► 57 decays. These backgrounds are suppressed by applying the selection requirements described 
below. These requirements have been optimized separately for each signal mode, using simulated 
signal and background event samples and the method described in [9], for maximum statistical 
sensitivity 2 assuming a branching fraction of 1.0(0.5) x 10~ 6 for the charged (neutral) mode. 

The photon from the signal B decay is identified as a localized energy deposit (cluster) in 
the calorimeter with energy 1.5 < E* < 3.5 GeV in the center-of-mass (CM) frame. The energy 
deposit must not be associated with any reconstructed charged track, be well-isolated from other 
EMC clusters, and meet a number of further requirements designed to eliminate background from 
hadronic showers, small-angle photon pairs, and charged particles [10]. 

We veto those photons that can be associated with another detected photon to form a it or r) 
candidate using the likelihood ratio 



r(M(^'),E 7/ \i) ; _ 



P(M( 7 7'), Ey |signal) + P(M( 7 y), E y \i) ' 



ttV (3) 



In this definition, V is the probability density function defined in terms of the invariant mass of 
the photon pair, 7^(77'), and the energy of 7' in the laboratory frame, Ey, as determined from 
simulated signal and background events. To consider photons coming from the decays of 7T° and 77 
that have converted to e + e~ pairs, we combine the high-energy photon candidate with any e + e~ 
pair in the event with an invariant mass m e + e - < 50 MeV/c 2 , and reject the photon if the total 
invariant mass satisfies either 100 < m^ e + e - < 160MeV/c 2 or 500 < m^ e + e - < 590MeV/c 2 . 

Charged pion candidates are selected from well-reconstructed tracks with a minimum momen- 
tum transverse to the beam direction of 100 MeV/c. In order to reduce backgrounds from charged 
kaons produced in b — > 57 processes, a 7r =l= selection algorithm [5] is applied, combining DIRC 
information with the energy loss measured in the tracking system. 

2 Here, the figure of merit is S/ \/{S + B), where S and B are the rates for signal and backgrounds respectively. 
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Photon candidates identified in the EMC with energy greater than 50MeV are combined into 
pairs to form 7r° candidates. For B° — > (B + — > /) + 7) decays, the invariant mass of the pair is 
required to satisfy 122 < to 77 < 150 MeV/c 2 (117 < to 77 < 148 MeV/c 2 ). We also require the 
cosine of the opening angle between the daughter photons in the laboratory frame be greater than 
0.413 and 0.789 for B° — ► loj and B + — > p + 7 respectively. 

The identified pions are combined into vector meson candidates by requiring 633 < m n + n - < 
957 MeV/c 2 , 636 < nv^o < 932 MeV/c 2 , and 764 < m^+^o < 795 MeV/c 2 for p°, p+, and cj 
respectively. The charged pion pair must originate from a common vertex, which is required to be 
consistent with the interaction region to suppress K® decays. 

The photon and p/uo candidates are combined to form the S-meson candidates. We define 
AE = E* B — E^ cam , where Eg is the CM energy of the S-meson candidate and ££ eam is the CM 

beam energy. We also define the beam-energy-substituted mass toes = y^beam ~~ Pi? 2 ' wnere Pjj 
is the CM momentum of the B candidate. Signal events are expected to have a AE distribution 
centered at zero with a resolution of about 50MeV, and a toes distribution centered at the mass 
of the B meson, tub, with a resolution of 3 MeV/c 2 . We consider candidates in the ranges —0.3 < 
AE < 0.3 GeV and toes > 5.22 GeV/c 2 to incorporate sidebands that allow the combinatorial 
background yields to be extracted from a fit to the data. 

To suppress B — > p(ir /r]) and B — > uj(tt°/t]) events, we calculate the vector meson helicity 
angle, 9h, defined as the angle between the n~ track (normal to the u decay plane) and the B 
momentum vector in the p (u>) rest frame. We require | cos 9h\ < 0.75. 

Contributions from continuum background processes are reduced by considering only events 
for which the ratio R2 of second-to-zeroth order Fox- Wolfram moments [11] is less than 0.7. In 
addition, several variables that distinguish between signal and continuum events are combined in a 
neural network. The quantity R' 2 , which is R2 in the frame recoiling against the photon momentum, 
is used to reject ISR events. To discriminate between the jet-like continuum background and the 
more spherically-symmetric signal events, we compute the angle between the photon and the thrust 
axis of the rest of the event (ROE) in the CM frame. The ROE is defined by all the charged tracks 
and neutral energy deposits in the calorimeter that are not used to reconstruct the B candidate. 
We also calculate the moments Lj = J2jPj ■ I cos Qj IV Pj: where p*j and 9j are the momentum 
and angle with respect to an axis, respectively, for each particle j in the ROE. We use L\, L 2 , and 
L3 with respect to the thrust axis of the ROE, as well as with respect to the photon direction. In 
addition, we calculate the l?-meson production angle 6* B with respect to the beam axis in the CM 
frame. Differences in lepton and kaon production between background and B decays are exploited 
by including flavor-tagging variables [12] as well as the maximum CM momentum and number of 
K ± and K® in the ROE. The significance of the separation along the beam axis of the S-meson 
candidate and ROE vertices is included as well. To reject events for which this quantity is poorly 
reconstructed, the separation along the beam axis and the associated uncertainty are required to 
be less than 4 mm and 0.4 mm, respectively. 

We train the neural network separately for each signal mode and select B + — > p + j, B° — ► p°j, 
and B° — > loj candidates with a requirement on the the neural-network output that retains 63%, 
74%, and 71% of the signal events respectively. For these cuts, we determine the continuum 
background efficiencies using a data sample of 27. 2 fb' 1 taken 40MeV below the T(4S) resonance 
as 3.0%, 5.3% and 6.7% for the three signal modes respectively. 

The expected average candidate multiplicity in the selected signal events is 1.01 for B° — > p°7 
and 1.07 for B + — > p + j and B° — ► lo^; in events with multiple candidates the one with the 
reconstructed vector meson mass closest to the nominal mass is retained. 
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Applying all the selection criteria described above, we find efficiencies of 11.6% for B + — > /? + 7, 
14.5% for B° -► p° 7 , and 8.1% for -► oj 7 . 

4 MAXIMUM LIKELIHOOD FIT 

The signal content of the data is determined by means of a multi-dimensional unbinned maximum 
likelihood fit, which is constructed individually for each of the three signal decay modes. All fits 
use AE, rriES , cos Oh, and the neural-network output AW, after transforming it according to 



V c 3 



constant (4) 



in order to facilitate the parameterization of the probability density function (PDF) used in the fit. 
For B° — > uj, the cosine of the Dalitz angle Od [5] is added as a fifth observable. 

In addition to signal and continuum background processes, we consider several sources of back- 
ground from B decays, which in the fit are combined in different ways depending on the signal mode 
under study. In the B° — > tvy fit, all B backgrounds are combined into a single component, while 
for the B° — > p°7 analysis B + — ► K* + j, B° — > -fT*°7, and other i? background processes are treated 
separately. The B + — > p + 7 fit uses four different categories of 5 backgrounds: 5 + — > i^T* + 7 with 
— > K + 7T , other 5 — ► if* 7 decays, B — > X s 7 processes (excluding 5 — ► if*7), and remaining 
.B backgrounds. 

In studies of simulated signal and background event samples, the correlations among the observ- 
ables are found to be small. We therefore assume that the PDF V(xj \ cti) for each of the A^yp event 
hypotheses is the product of individual PDFs for the fit observables Xj given the set of parameters 
Si. The likelihood function for signal mode k (= p + 7, p°7, toj) is defined as 



v hyp 



7V h; 

C k = exp ( - n i 

i=i 



' N k /N hyp 
3=1 \ i=l 



(5) 



where rii is the yield of each hypothesis and is the number of candidate events observed in data. 

The functional form of each PDF is determined from a one-dimensional fit to a dedicated sample 
of simulated events. The AE distribution is corrected for the observed difference between data and 
simulated samples of B — > K*^/ decays. All continuum background PDF parameters float freely 
in the fits while the shapes of the signal and B background distributions are fixed. For B° — > ujj, 
the B background yield floats freely in the fit. In the B + — > p + 7 analysis, the B + — > K* + ~f 
(K* + — ► i^ + 7r°) contribution and the ratio of the other three -B background yields are determined 
from simulated events, as are the relative contributions from the three B background components 
in the B° -> p°j fit. 

For the signal, the tties spectra are described by Crystal Ball functions [13], the angular dis- 
tributions are modeled by second-order polynomials, and the distributions of AE and MM are 
parametrized as asymmetric, variable- width Gaussians 



f(x) = exp 



-{x-pf 



2a L,R + OClM x ~ 



(6) 



where [i is the peak position of the distribution, ol,r are the width left and right of the peak, and 
(Xl,r are a measure of the tail on the left and right side of the peak respectively. 
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The function © also describes the continuum background MM shape; the remaining continuum 
spectra are modeled by ARGUS functions [14] (m-Es) ° r second- and fourth-order polynomials (AE, 
cos 8h, and cos#d). Various functional forms are used to describe the different B background 
components. 

In order to measure the combined branching fraction B[B — > (p/u)j], we also perform a si- 
multaneous fit to the three decay-mode specific data sets for the effective signal yield n e fr, which 
is related to signal yields and reconstruction efficiencies 3 obtained from the individual fits via 
n (£+ ^ p + 7 ) = neff . l e ( B + ^ p + 7 ) and n ( B ° ^ (pP/ w ) 7 ) = i^n eff • e(B° -» (p%j) 7 ). 

Figures |H [3 and ^ show the projections of the fit results for B + — ► p + 7, 5° — > /0°7, and 
-B — > a;7 respectively compared to the data; for each plot the signal fraction is enhanced by 
selections on the other fit variables. The resulting signal yields are given in Table[21 The significance 
is computed as \J2A log C, where A log C is the log-likelihood difference between the best fit and a 
fit to the null-signal hypothesis; only statistical uncertainties are included here. 

5 SYSTEMATIC UNCERTAINTIES 

Tabled gives an overview of the contributions to the systematic uncertainties. These are associated 
with the signal reconstruction efficiency, the modeling of BB backgrounds, and the choice of fixed 
parameters of the fit PDFs. The latter two contribute to the uncertainties on the signal yields. A 
small uncertainty on the overall normalization is associated with the imperfect knowledge of the 
total number of BB pairs in the underlying data sample. 

Table 1: Fractional systematic errors (in %) of the measured branching fractions. 



Source of error 


B + - 


> p+7 


B° -> A 


B° 


Tracking efficiency 


1.0 


% 


2.0% 


2.0% 


Charged-particle identification 


2.G 


i% 


4.0% 


2.0% 


Photon selection 


1.9 


>% 


2.6% 


1.7% 


7T° reconstruction 


3.0 


'% 




3.0% 


7T° and r\ veto 


2.8 


% 


2.8% 


2.8% 


MM efficiency 


5.0 


% 


3.5% 


3.5% 


AE shape from K*j 


3.1% 


2.4% 


1.9% 


MM shape 


0.2% 


3.9% 


4.7% 


B background normalization 


3.0 


i% 


4.0% 




B counting 


1.1% 


1.1% 


1.1% 


Combined 


8.4 


% 


9.2% 


8.2% 



The signal efficiency systematic error includes uncertainties from tracking, charged-particle 
3 The efficiencies include the daughter branching fractions. 
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Table 2: The signal yield (n S i S ), statistical significance in standard deviations (a), efficiency (e), 
and branching fraction (£>) central value for each mode. The errors on (n S i g ) are statistical only, 
while for the branching fraction the first errors are statistical and the second systematic. All results 
are preliminary. 



Mode n sig Significance e(%) £(1(T 6 ) 

B+ _> p + 1 42.4±i£S 4.1cj 11.6 1.06±gjf ± 0.09 

B° -> p° 7 38.7i^° 8 6 5.2a 14.5 0.77±g;f^ ± 0.07 

B°^coj H-Olsle 2.3<T 8.1 0.39^;^ ± 0.03 (< 0.84 at 90% C.L.) 



identification, 7/71" reconstruction, photon selection and the neural network selection that are 
determined from suitable independent data control samples. 

To estimate the uncertainty related to the extraction of the signal PDFs from MC distributions, 
we vary the parameters within their errors. The uncertainty related to the choice of a specific 
functional form for the shape of the MM distributions is evaluated by using a binned histogram as 
an alternative PDF. All relative and absolute normalizations of B background components which 
were fixed in the fit are varied by 50%. For all these variations, the corresponding change in the 
fitted signal yield is taken as a systematic uncertainty. 

6 RESULTS 

The branching fractions are calculated from the fitted signal yields assuming B(T(4S) — ► B°B°) = 
B(T(45) — > B + B~) = 0.5. For B° — > ujj, we also compute the corresponding 90% confidence 
level (C.L.) upper limit using a Bayesian technique. The signal yield upper limit n/ is determined 
such that Jq 1 Ldnj / °° Cdn = 0.90, assuming a flat prior. The systematic uncertainty is included 
by increasing ni and decreasing the detection efficiency by their respective errors. The results are 
listed in Tabled 

The simultaneous fit finds an effective signal yield n e ff = 702^^^ with a corresponding statistical 
significance of 6.3 a. This translates into a combined branching fraction 

B[B -> (p/u)i\ = (1.01 ± 0.21 ± 0.08) x 10~ 6 . (7) 

We also measure the ratio T(B + p + j)/[2T(B° p°>y)} - 1 = -0.36 ± 0.27 in order to test the 
hypothesis of isospin symmetry. The result is in agreement with the theoretical expectation [2]. 
Using the measured value of B{B — > K*j) [15], we calculate 

B[B -» (p/u)j]/B(B -> K*~f) = 0.024 ± 0.005. (8) 

This result is used to determine the ratio of CKM elements |T4d/^ts| by means of Equation |2J). 
Following [16], we choose the values l/( = 1.17 ± 0.09, and AR = 0.1 ± 0.1. We find 

It/ /t/ I _ n 171 +0.018+0. 017 / Q \ 
Ind/^tal - U-171_ .o21-0.014> i 9 ) 
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where the first error is experimental and the second is theoretical. This is consistent with the 
current world average of \V td /V ts \ = 0.201±g;ggf [17]. 

Using the measured value of B(B° — ► K * 7) [15], we also calculate 

2 x B(B° -» p°~f)/B(B° -> iT° 7 ) = 0.038±Hio- ( 10 ) 

By only using these two neutral decay modes, the theoretical interpretation of |VtdA^s| is simplified 
since the VF-annihilation processes present in the B + — > p + 7 channel are avoided. Analogous to 
Equation (j2j), taking the same values for 1/Q and AR as above, this result is used to obtain 

\Vtd/V ts \ p o/ K *o = 0.216lo:o3iiaoi8> ( n ) 
where the first error is experimental and the second is theoretical. 



7 SUMMARY 

In conclusion, we observe the exclusive b — ► dj transitions B + — > p + j and B° — » p°7 and measure 
the branching fractions £(£+ -> p+ 7 ) = (1.06±g;jjf ± 0.09) x 10~ 6 and B(B° -» p°7) = (0.77lg;S ± 
0.07) x 10 -6 , where the first error is statistical and the second is systematic. We set an improved 
90% C.L. upper limit on the B° — ► W7 branching fraction of B(B° — ► CJ7) < 0.84 x 10 -6 . Assuming 
isospin relations between the three branching fractions, we measure the combined branching fraction 
B[B — ► {p/u))^\ = (1.01 ± 0.21 ± 0.08) x 10~ 6 . This result translates into a measurement of the 
CKM matrix element ratio |Vt<f/^4s| = 0-17llo'o2i( ex P-)-o oi4(^ eor ')- ^ n addition, we measure the 
isospin asymmetry Y(B + — > p + 7)/[2r(i?° — > p ^/)] — 1 = —0.36 ±0.27. All these preliminary results 
are consistent within errors with the SM predictions. 
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Figure 2: Projections of the fits to the B + — ► p + ^f sample in the discriminating variables AE (upper 
left), tties (upper right), J\fJ\f (lower left), and cos Oh (lower right). The points are data, the solid 
line is the total PDF and the dark dashed (light dot-dashed) line is the background (signal) only 
PDF. The selections applied, unless the variable is projected, are: —0.15 < AE < 0.05 GeV, 
5.275 < m ES < 5.285 GeV/c 2 , and AfAf > 0.0. 



16 




A E (GeV) 



m ES (GeV/c') 




-2 2 4 6 8 

p° transformed NN output 



o 30 

£ 25 
c 

CD 

> 20 

LU 

15 
10 
5 




BaBar 

preliminary 




-0.6 -0.4 -0.2 



0.2 0.4 0.6 
C0S Helicity 



Figure 3: Projections of the fits to the B — ► p°j sample in the discriminating variables AE (upper 
left), niES (upper right), NN (lower left), and cosOh (lower right). The points are data, the solid 
line is the total PDF and the dark dashed (light dot-dashed) line is the background (signal) only 
PDF. The selections applied, unless the variable is projected, are: —0.15 < AE < 0.05 GeV, 
5.275 < m E s < 5.285 GeV/c 2 , and NN > 0.0. 
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Figure 4: Projections of the fits to the B — * u>j sample in the discriminating variables AE 
(upper left), tues (upper right), MM (middle left), cosOh (middle right), and cos9d (bottom). 
The points are data, the solid line is the total PDF and the dark dashed (light dot-dashed) line 
is the background (signal) only PDF. The selections applied, unless the variable is projected, are: 
-0.15 < AE < 0.05 GeV, 5.275 < m E s < 5.285 GeV/c 2 , and MM > 0.0. 
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